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1. Introduction 

Several groups are investing considerable effort to numerically calculate various 
hadronic decay constants and the mass spectrum with Wilson fermions [Q]. One 
of the systematic errors, that due to the finite lattice spacing, can in principle be 
diminished by taking an improved action a la Symanzik . In the weak coupling ex- 
pansion the pure Yang-Mills theory has no 0(a) cutoff effects in the infinite volume 
or with a finite volume with periodic boundary conditions, however the fermionic 
part of the Wilson action introduces such terms. To eliminate these effects all inde- 
pendent dimension five operators are added to the action. The procedure drastically 
simplifies by demanding that only the on-shell physical quantities of the theory are 
improved Then only the operator 'tl){x)a fj_yPfj_y{x)%l){x) needs to be added to the 
action where Pfj^p is the term with four plaquettes touching at point x in the form 
of a 'clover leaf , an equivalent of F^j^y 0]. The fermionic part of the action reads 

= a^Xl [(p{x)[{r - -i^)U^{x)^{x + iia) + {r + -i^)U^\x - ^a)ii{x - ^a)] 

(1) 

with 

^^.(^) = IKA^) + f/4.(x) + U^,_,ix) + - ifi^ u)] (2) 

where Ul^,^{x) is the untraced plaquette variable. The coefficient c^Qq) ought to 
be fixed perturbatively to all orders of Qq. This coefficient happens to be quite 
important for the hadronic decay amplitude and for the hadron spectrum when the 
short distance effects play a role e.g. for the hyperfine splitting in charmonium. At 
tree level, the improvement condition fixes this to be 'r' the Wilson parameter |Q. 
The one-loop coefficient was first calculated by WohlertQ and also there is a mean 
field estimate given by Lepage and Mackenzie [^. In this article we independently 
check this coefficient at the one loop level and give an approximate value in the two 
loop level. 

2. General Strategy 

To fix the coefficient c((7o) to all orders of perturbation theory a suitable on-shell 
quantity needs to be defined. The most obvious of such quantities which instantly 
comes to mind is to look for the pole of the fermion two point function. This however 

^ The action given in Ref. |^ is not exactly the clover action now used in numerical simulations. 
We found terms missing from the fermion-three gluon vertex function and errors in the expression 
of Clebsch-Gordon coefficients. This created a suspicion of a drastic change of the result since it 
involves the most dominant tadpole graph. We made an independent check however and found 
that the number quoted is correct and these errors are probably just misprints. 
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unfortunately cannot be related to c((7q) due to the Slavnov- Taylor identity 

-z^s(p) = r^(p,p,o), (3) 

and the fact that the improved part of the vertex function T^{p, p, k) (oc J2u ^^f^u sinky) 
vanishes for = 0. This clearly shows that to all orders of perturbation c{gQ) is not 
sensitive to the on-shell fermion two point function. Another quantity of interest 

is < 'il){y)Pe^°^ y^^^^^'^^ 'tp{x) >, however does not yield a tree level improvement 
condition. We were so far unable to find a suitable quantity which needs an im- 
provement at tree level in the infinite volume. However, one possibility is to assume 
that our lattice has finite extent L in xi and X2 directions with twisted periodic 
boundary conditions. Due to these boundary conditions, quarks and gluons get 
mass by a Kaluza-Klein type of mechanism and some of the eigenvalue modes of the 
transfer matrix remain stable for small coupling in this twisted world. These states 
are created from the vacuum by gauge invariant operators such as a Wilson loop 
winding around the torus. The spectrum of the S-matrix for these gauge invariant 
modes in the LSZ scheme can be studied unambiguously as has been already done 
for the pure Yang-Mills theory to fix two necessary coefficients 0. A more detailed 
discussion of this procedure is given in Ref . . 

The twisted periodic boundary condition for the gauge field reads, 

U{x + L^, /i) = n,U{x, fi)n,-\ z/ = 1, 2, (4) 

where Qi, are constant SU(N) matrices with 

Q^Q^ = e^n2^i- (5) 

To use the twisted boundary condition for fermions one introduces an extra internal 
degree of freedom "smell" and a fermion belongs to the Ns x Nc representation. 
The twisted antiperiodic boundary condition on the torus reads 

ip{x + Ly) =Qy'il){x)Qy~^e'^ z/=l,2. (6) 

Thus the transverse momenta in these directions take discrete values as 

p± = {2ni + 1, 2n2 + l)m, rii G Z, (7) 

where m = This offers a mass gap even though we have a zero bare mass 

fermion to start with. The Fourier decomposition of the gauge field reads 

A^{x) = -^Y^k^j exp{ikx) Tk exp ^i^(A;) (8) 

where k±_ = 2mny and plays the role of SU(N) group generator. Since is 
traceless, there exists no zero mode. The gluon propagator reads 

D^^uik) = -\ Xk Z{k, k)[5^,P + {a- l)kM-^^, (9) 
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where 



Xk= z/ A;x = 0, (modN) 
= 1 otherwise, 

(10) 

a is the gauge fixing parameter and = 2 sin^. Also 

Z(A:,A;') =^^{<'^''='>-('=''=')) (n) 

where z = ,< k, k' >= riirii + 72,2712' + {rii + 712) {ni + ^2') and {k, k') = 71172,2' — 
n2ni . Following Ref. [Q we study the spectrum of the LSZ scattering process 
(c.f. Fig.l). The on-shell momenta for these fermions are pi = {iE,m,m,im) , 
Pi = {iE,m, —m, —im) , p2 = {iE,m, —m, ~im) and ^2' = {iE,m,m,im). This 
kinematical choice of momenta simplifies the calculation drastically since the ex- 
change diagram Fig. lb does not contribute and also the fermion wave function 
renormalization to 0(a) does not contribute to the S-matrix elements. Following 
Ref . , we look for the residue of the pole of the scattering amplitude 

S = T^ZUk)D^,{k)T,', (12) 

where = U{pi) T^{pi,pi', k) U{pi) and Zcik) is the gluon wavefunction renor- 
malization. Also Tj can be obtained by replacing pi P2 , Pi' P2 and k —k 
in T^. The residue of the pole of this S-matrix element; 

\ = z{-k,p^')T^T;. (13) 

We assume here a perturbative expansion of this residue A = YJ^i f/o^* A* and also 
forT^ = f?o''~^' V and c(f72) = Cq + glc^ +9^02 + . . .. 

To illustrate the case in the tree level 

= T;il) + T'^ia) = zik,p)U^ip)[~^,~'-{co-r)ip+p'), + Oia')]U,{p'). (14) 

We look for the residue in a particular polarization of the gluon say fi = 1 and 
in the fixed helicity state of the fermion say a = (3 = 1. For the zero bare mass 
fermion [/(l?) = (sinh i?)~^.(7o sinh i? + i^iPi) ?7^, where [/^ = 5i^p. This gives 

T°(l) = 2iz, T^{a) = -2izm{cQ-r) 
r°(l)'(l) = -2i, T°(a)' = 2im{co-r) 

(15) 

and 

A°(a) = (T^il) n\a) + ^{0)^' {I)) 
= — 8zm{co — r). 

(16) 
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To demand the tree level improvement we set A°(a) = and thus get the condition 
Co = r i. 

To all orders of perturbation theory Ti(a) gets contributions from the fermion 
wave function renormalization, vertex function renormalization and also a term com- 
ing from the naive expansion of c(5'q). So 

ri(a) = +rf^ +r,^^. (i7) 

However for our kinematical choice of momenta, it can be proved that 

T°(l) ^'(a) +Tf^(a)Tf(l) =0. (18) 

The contribution of the gluon wave function renormalization Zdk) does not con- 
tribute to A for the 0{a) improvement. So it remains to calculate only T^^{a) and 
T^^'{a) to all orders of perturbation. The three point function U (p) r^(p, p', k) U {p') 
is expressed as Yl,i Oi Bi where Oj's are sixteen bilinear invariant Dirac basis 
matrices and Bi 's are their coefficients. It is sufficient to look for the residue in the 
particular channel which we have chosen, here oc ai2 for every diagram. 

3. 0{g^) Improvement 

To one loop order there are six diagrams (c.f.Fig.2) contributing to Ti(a). Due 
to the twisted periodic boundary conditions all these loop integrals depend on L 
and N. Then the contributions are compared with their asymptotic expansion 

^~E.=o("^ 1^ m)(m)^+'=' (19) 

where the lattice spacing is set to unity, = yn ^ degree of divergence 

of the graph. The method of evaluation of these graphs is exactly the same as in 
Ref . . The analytical expressions for the vertex functions and loop integrands are 
quite lengthy and beyond the scope of this paper . So we present here only the 
results. 

The imaginary part and the coefficient of the log m for daigram 2.b exactly 
cancel with that of 2.c. This is also true for the diagram 2.d where this part cancels 
with the sum of the contributions coming from 2.e and 2.f.0 The small contributions 
coming from diagram 2.b-f are not checked by us, however we use here the results 
of Wohlert The contribution of each diagram is given in Table I. for N = 2 and 
= 3 in the Feynman gauge. Needless to state that the total contributions of these 
graphs to the residue of the S-matrix is independent of the choice of the gauge. 

Table 1. Values of the one-loop graphs for N = 2 and = 3. 

^ This can be analytically checked in the following way. We take a very small generic external 
momenta e for these diagrams and then make a Taylor's expansion in e of these loops to get 
the same tensor structure of all these loops by using symmetry of the integration of the internal 
momenta. (In the Feynman gauge expressions are considerably simpler.) Then we look for the 
logarithmic divergence term for each loop. 
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Fie;. 


N=2 


N=3 


2. a 


0.34800 


0.594789 


2.b 


-0.01630 


-0.03268 


2.C 


0.01512 


0.03087 


2.d 


0.03864 


0.05787 


2.e 


-0.03903 


- 0.05953 


2.f 


-0.03903 


- 0.05953 


Total 


0.308 


0.53179 



Here also one observes that the contribution of 2.b and 2.c nearly cancel with each 
other so also of 2.d nearly cancel with 2.e .The only leftovers are the contributions 
of 2. a which is the most dominant one and the contribution of 2.f is nearly 10% of 
2. a. This clearly shows the evidence of the tadpole dominance ^j. 

So the residue 

X\a) = 4m [0.308 - 2ci], (20) 

for iV = 2 and 

X\a) = Am [0.5318 - 2ci] (21) 
for = 3. This gives Ci = 0.154 for = 2 and Ci = .2659 for = 3. 

4- 0{g^) Improvement 

The numerical simulations for the hadron spectrum or decay amplitudes are 
performed for ^ 6. i.e. ~ 1. So it is worthwhile also to fix the coefficient C2 
. From our observation in the previous section and also the earlier investigation of 
perturbative computation of lattice graphs |0| there is an indication that the main 
contribution comes from the tadpole graph. In our case at one loop level all graphs 
nearly cancel with each other except the tadpole one and also all are less than 10% 
of the latter. Thus it seems legitimate to take only the perturbative correction from 
tadpoles as a first step towards a full and much lengthier calculation of the two loop 
improvement. The dominant two loop graphs are given in Fig. 3. and the results 
are presented in Table II. Here again one also observes the tadpole dominance. 
The main contribution comes from Fig.3.a and the gauge invariant part of Fig. 3c. 
The contribution of the latter which is like plaquette-plaquette correlation function 
Di_iu,^iu, gives 0.1489, the largest among all these graphs. All other contributions are 
less than 10% of these two quantities. 
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Table II. Values of the two-loop graphs for N = 2 and = 3. 



Fig. 


N=2 


N=3 


3. a 


0.03403 


0.10465 


3.b 


-0.0259 


0.00648 


3.(c+d+e+f+g) 


0.05346 


0.16174 


3.(h+i) 


-Uf 0.01 


-Uf 0.02 



This gives ca = 0.0207 for A = 2 and ca = 0.1164 for =3 and nj = 2. 

5. Conclusion 

A meanfield type of estimate of this coefficient is given by]^ 

c{9o) = {^<trUp>)-i. (22) 
The perturbative expansion of this plaquette expectation value reads P 

^ < trUp >= 1 - uigl - u^gl + . . . (23) 

where 

(Ar2-l) (Ar2-l) 2A2-3 , , , 

with Ki = — 0.0042. This gives C\ = .25 and C2 = .098 whereas we get 
Ci = .266 and C2 = .1164 with two flavors. It is quite gratifying to observe that 
by two different methods of calculation one gets a fair agreement of the estimate 
of these coefficients, which can now be used in numerical simulations with more 
confidence. 
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Figure captions 

Fig.l. The on-shell fermion fermion scattering 

Fig. 2. The one- loop graphs 

Fig. 3. The most dominant two-loop graphs 
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